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ABSTRACT 

Using far-infrared imaging from the "Herschel Lensing Survey" , we derive dust properties of 
spectroscopically-confirmed cluster member galaxies within two massive systems at z^0.3: the 
merging Bullet Cluster and the more relaxed MS2137. 3-2353. Most star-forming cluster sources 
(~90%) have characteristic dust temperatures similar to local field galaxies of comparable infrared 
(IR) luminosity (T dust ~ 30 K). Several sub-LIRG (L m < 10 11 L ) Bullet Cluster members 
are much warmer (T^ust > 37 K) with far-infrared spectral energy distribution (SED) shapes 
resembling LIRG-type local templates. X-ray and mid-infrared data suggest that obscured active 
galactic nuclei do not contribute significantly to the infrared flux of these "warm dust" galaxies. 
Sources of comparable IR-luminosity and dust temperature are not observed in the relaxed cluster 
MS2137, although the significance is too low to speculate on an origin involving recent cluster 
merging. "Warm dust" galaxies are, however, statistically rarer in field samples (> 3 a), indicating 
that the responsible mechanism may relate to the dense environment. The spatial distribution 
of these sources is similar to the whole far-infrared bright population, i.e. preferentially located 
in the cluster periphery, although the galaxy hosts tend towards lower stellar masses (M* < 
10 10 M Q ). We propose dust stripping and heating processes which could be responsible for the 
unusually warm characteristic dust temperatures. A normal star-forming galaxy would need 30- 
50% of its dust removed (preferentially stripped from the outer reaches, where dust is typically 
cooler) to recover a SED similar to a "warm dust" galaxy. These progenitors would not require a 
higher IR-luminosity or dust mass than the currently observed normal star-forming population. 

Subject headings: Galaxies: clusters: general - galaxies: star formation - infrared: galaxies 
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1. Introduction 

High energy photons (i.e. ultraviolet light) heat 
dust grains and are re-emitted at longer wave- 
lengths, such that far-infrared (FIR) observations 
are an effective probe of dust properties. In galax- 
ies lacking an IR-bright active galactic nucleus 
(AGN), the dominant heat source are young stel- 
lar populations, and hence FIR flux also traces 
star formation (SF). The spectral energy distri- 
bution (SED) of a dust cloud can be modeled by 
a composite of many modified blackbodies (grey- 
bodies) , as the radiative temperature of each grain 
is a function of size and composition. With coarse 
photometric data, the dust component is often 
treated as a single greybody at a fixed characteris- 
tic temperature and emissivity, encapsulating the 
typical dust p roperties of the galaxy (see discus- 
sions m e.g. Ifflldebrandl Il983t iNetterfield et al 
20091) . 



Observing in the near and mid-infrared, Spitzer 
glimpses the Wien's law (blue-ward) side of the 
dust component. Generally, MIPS 70 and 160 
/mi data are too shallow to detect cluster galaxies 
other than the nearest a nd brightest (e.g. Sha pley 
Supercluster at z<0.05; lHaines et al.ll2011al) . In 
contrast, the superb sensitivity of MIPS at 24 /itm 
has been the de rigueur probe of obscured star for- 
mation. Using a small number of local luminous 
infrared gal a xies ( LIRGs) as a template library, 
Rieke et al. (2009) derived a simple, yet power- 



ful, formula to extrapolate obscured star forma- 
tion rate (SFR) directly from 24 /im flux alone. 

Cluster member galaxies are typically quies- 
cent early-types with no significant star fo rmation, 
as de picted by the m orphology-density dDresslerl 



So while SF does concentrate in pockets on the 
periphery of clusters, the star formation rate of 
star-forming cluster galaxies is no different from 
i dentically selected galaxies at lower densities 
(jGeach et al.l 120111) . Similarly, elevated number 
counts, but no enhancement in individual SFRs, 
are reported in the most extre me environments, 
such a s massive cluster merge rs ( De Propris et al 



2008; 



Chung et al 



2007t Uohnston-Hollitt et al 

I2010I ). Once star formation in a cluster member 
ceases, most likely due to gas exhaust ion before 



comp leting the first in-fall pass (e.g. iTreu et al 



2003), further gas cooling is prevented by the hot- 



ter intracluster medium of the dense core, giving 
rise to the observed density relations. An obvious 
exception are the brightest cluster galaxies at the 
center of relaxed 'cool-core' clusters, which exhibit 
star formation rates up to ~ 100 M w yr~* fuel ed 
by inflowing cool cluster gas (jRawle et al.ll2012i ). 



The Herschel Space Observatory (jPilbratt et al 



2010h directly constrains the full FIR dust com- 
ponent for a large sample of sources. A small 
number of early studies analyzed the FIR sig- 
natures of cluster galaxies at z^0.2-0.3. Wide- 
field analysis recovered a peak in star formation 
rate at 



-2-2.5 Mpc; ISmith et al.1l2010l 
Pereiraet all l2010h . However, the real power of 



Herschel is in the constraint of du st character- 
istics beyo nd merely IR lumin osity. iRawle et al 



1985c iPoggianti et al 



(|2010h and lPereira et al.l (|2010h reported an unex- 
pected cluster population displaying a FIR SED 
shape unlike local templates. Specifically the 
sources exhibited excess flux blue- ward of the dust 
peak, but were not associa ted with infrared-bright 
AGN. Rawle et al.l (|2010l ). an analysis of science 
demonstration phase (SDP) Herschel observations 



1999) relations. Ultravio- 



1980) and SF-densitv ^Dressier. Thompson fc Schechtmanl of the massive, merging Bullet Cluster found that 

^40% of star-forming cluster galaxies have ele- 
vated Sioo/SW compared to local templates. Sim- 
ilar sources were not identified in field surveys at 
com parable redshift, or in high-redsh ift samples 
(e.g. lElbaz et alj|20ld; iRex et al.|[201oh . 

In the current study, far-infrared observations 
from the Open Ti me Key Program, t he ll Herschel 
Lensing Survey" ( Egami et al.ll201o( ) are used to 
explore the dust properties of galaxies in two mas- 
sive clusters at z^ 0.3. The paper has the following 
structure: Section [2] introduces the target clusters 
and data; Section[3]describes photometric analysis 
and SED fitting; Section |4] explores the IR lumi- 



let observations do, however, reveal low-level re- 
cent star formation in s ome optical-red-seq uence 
sources (|Yi et al.l 120051 IRawle et all 120081 ). In- 
frared imaging indicates that galaxies (in small 
groups) accrete onto a cluster along filamentary 
structures, where gravitational interactions be- 
tween the galaxies, rather than clus ter-potential 



processes, stimulate star f o rmation (IFadda et al 



2008; iKovama et ail 120081 lHaines et all l2011bf ). 



space observatory with science instruments provided by 
European-led Principal Investigator consortia and with im- 
portant participation from NASA. 
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nosity, dust properties and physical interpretation; 
Section [5] summarizes the primary conclusions. An 
in-depth examination of the SED fits is detailed in 
the Appendices. 

2. Observations 



to ram pressure (Clowe. Gonzalez fc Markevitch 



Th e " H erschel Lensing Survey" (HLS: lEgami et al 
2010) is designed to exploit the lensing effect 
of massive clusters, pushing beyond the nomi- 
nal confusion limit to observe intrinsical l y faint , 
high-redshift source s, e.g. iRex et ah (l2010h . 
Co mbes et al. (2012). To this end, H LS comprises 
deep PACS ([Poelitsch et al.1 l2010h and SPIRE 
(jGriffin et al-l boiOl) imaging (100-500 /mi) of 44 
massive cluster cores, 0.2 < z < 0.5. 

In this analysis, we concentrate on galaxies 
within two massive systems: the Bullet Clus- 
ter and MS2137.3-2353. These clusters have 
substantial multi-wavelength photometric cov- 
erage, alongside 100s of spectroscopic redshift 
measurements. Although photometric redshifts 
derived from optical-FIR photometry are suffi- 
ciently accurate to identify e.g. the high red- 
shift population ( Perez-Gonzalez et al. 2010l ). iso- 
lation of cluster members requires precise spec- 
troscopic redshifts. Furthermore, redshift and 
dust temperature are degenerate when examin- 
ing the FIR SED. By placing a 30 K greybody at 
z~0.3, and then assuming various redshifts within 
the typical photomet ric redshift error (Az~0.1; 
Perez-Gonzalez et all ), we estimate the associated 
temperature uncertainty to be >10 K. 

In this section we detail relevant observations, 
while Figure [T] displays the spatial coverage of the 
infrared imaging. 

2.1. The Bullet Cluster (IE 0657-56) 

iRawle et all (|2010h presented a preliminary 
analysis of the HLS SDP data available for the 
Bullet Cluster (RA = 06:58:31.10, Dec = 
55:56:49.0; z = 0.296 ). The system is a recen t col- 
lision of two clusters (jMarkevitch et aDl2002 ), per- 
pend icular to the line of sight (jMarkevitch et al.l 
20041 ) . offering a unique laboratory for the study 
of SF within a dense dynamic environment. X- 
ray emission shows a supersonic bow shock pre- 
ceding the hot gas, while the weak lensing mass 
profile indicates that this X-ray bright compo- 
nent lags behind the sub-cluster galaxies due 



2004; IClowe et al.l 120061 ). Ram pressure from the 
merger event does not appear to impact the star 
formation rate of galaxies in the vicinity of the ho t 
X-ray gas (jChung et al.ll2009l; IRawle et al.ll2010l) . 

As described in the introduction, early Herschel 
cluster observations identified a surprising popula- 
tion of sources with an apparent excess at 100 /im. 
We obtained additional PACS coverage of the Bul- 
let Cluster to explore these sources further. 

2.1.1. Herschel/PACS observations 

The SDP PACS observations (pr esented in 
Egami et al. 20101 : Rawle et al. 2010l) consisted 
of two orthogonal scan maps, each comprising 18 
repetitions of thirteen parallel 4' scan legs. Total 
observing time was 4.4 hours, with 1500 sec/pixel 
integration time. PACS images are produced us- 
ing a modified versio n of the s tandard reduction 
package, HIPE v6.0 (|Otdl2010l) . High-pass filter- 
ing of the time-stream is necessary to remove the 
1/f noise drift of the PACS bolometers, and we 
implement a simple masking algorithm to avoid 
ringing of bright sources. Raw PACS data in- 
cludes a large number of usable frames beyond the 
science scan legs (turnaround data). We include 
all frames with speeds within 5% of the nominal 
scan speed. As the detector footprint is compara- 
ble in size to the final map, the additional frames 
increase signal-to-noise in the central regions, as 
well as overall spatial extent. Thus, SDP coverage 
of the central 8'x8' reached a 3<r depth of 3.1 and 
5.8 mJy in the 100 and 160 /im bands. 

The new, augmented PACS imaging is an ex- 
tension on the SDP data in three senses: (1) the 
100 and 160 /im footprints are approximately dou- 
bled (~8'xl5') by a new observation (identical 
configuration), to the south-west of the original 
pointing, chosen for the number of potentially in- 
teresting sources identified in the larger SPIRE 
images; (2) both pointings were observed in a sin- 
gle, deep 70 /tm observation, probing the SED be- 
tween existing 24 and 100 /tm data; (3) the 70 /tm 
pointing simultaneously observes 160 /tm, provid- 
ing extra depth in the longest PACS band. The 
70 /im map consists of two orthogonal scan maps 
of 21 repetitions each. The first has 28 parallel 
scan legs of 5.2' length, while the second has six- 
teen parallel 9' scan legs, with the pointing centers 
chosen to create a rectangular map (^15' x 8' with 
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Fig. 1. — 20'x20' fields aligned to the nominal cluster centers: Bullet Cluster (left) and MS2137 (right). Imag- 
ing footprints are shown as outlines: IRAC=grey, MIPS 24 /im=green, PACS=orange, SPIRE=magenta. 
Spectroscopically confirmed cluster members are indicated: crosses are undetected at A > 24 /im; unfilled 
symbols are detected at MIPS 24 /im but not by Herschel; filled symbols are detected by Herschel. Initiating 
the convention for plots in this paper, Bullet Cluster=red circles, MS2137=blue squares. The larger number 
of Bullet Cluster detections results from more substantial spectroscopy. 



turnaround data). Total observing time is 10.5 
hours, with ^2100 sec/pixel integration time. The 
final 70 /im map is actually larger than the two 100 
/im footprints, resulting in a significant number of 
sources near the edge of the PACS coverage with 
70 and 160 /im photometry only. 

The three final maps reach mean 3 a detection 
limits of 2.7, 2.8, 3.8 mJy at 70, 100, 160 /im re- 
spectively (beam sizes of 5.2, 7.7, 12"; detection 
limits calculated using the method described in 
IPerez-Gonzalez et al.l[2010h . These maps are all 
abov e the (3 a) confus ion limits of 0.07, 0.3, 2.3 
mJy (jBerta et alJl2010h . 



2.1.2. Herschel/ 'SPIRE observations 

The Bullet Cluster has not been re-observed by 
SPIRE since SDP, which used large map mode: 
20 repetitions each with a nominal length and 
height of 4'. The total observing time was 1.7 
hours, with 17 sec/pixel. The images were pro- 
duced via the standard reduction package HIPE 
v5.0, but include all turnaround data with speeds 
>0.5"/s to increase signal-to-noise in the outer 



regions. Final maps are complete to a cluster- 
centric radius >9', with 3er detection limits of 
8.2, 9.6 12.4 mJy at 250, 350, 50 /im respec- 
tively ( Perez-Gonzalez et al.l |2010( ). All three 
bands are confusion limit ed (3 <7 con f ~ 17, 19, 20 
mJy; iNguven et al. 2010l ). with beam sizes of 18, 
25, 36". 

2.1.3. Additional imaging 

As a consequence of the unique cluster char- 
acteristi cs, and a well known , lensed background 



source ( Gonzalez et al. 2009), there is a wealth 



of available ancillary data for the Bullet Clus- 
ter. The central 13'xl3' (including the entire 
PACS field) is covered by 5-band Spitzer imaging 
([Gonzalez et al.l l2009h from IRAC (3.6, 4.5, 5.8, 
8.0 /xm) and MIPS (24 /im), with 3<j detection 
limits of 2, 3, 5, 8 and 50 /iJy respectively. Also 
available are deep wide-field images i n BVR from 
Magellan/IMACS (|Clowe et al.ll2006l) . YJH from 
VLT/HAWK-I (PI: J.-G. Cub y) and X-ray from a 
0.5 Msec Chandra in tegration (iMarkevitch I l2006h . 
See the appendix of lEgami et al.l 2010l for further 
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details. 

2.1.4- Spectroscopic coverage 

Spectroscopic redshifts in the Bullet Cluster 
field originate from three campai gns: Magel- 
lan I MACS multi-slit (856 spectra; IChung et al. 
2010h . CTIO Hydra multi-fiber (202; D. Fadda et 
al. in preparation) and VLT FORS multi-slit (14; 
J. Richard, private communication). The merged 
catalog comprises 929 sources within the SPIRE 
field, including 371 cluster members (0.28 < z < 
0.31). 47 cluster galaxies are included in multi- 
ple spectroscopic datasets, with no membership 
disagreement. Figure [T] indicates the location of 
spectroscopic cluster members with respect to the 
infrared imaging coverage. 



sensitivities as in the Bullet Cluster. The entire 
SPIRE field is enclosed by deep BVRIZ Subaru 
Suprime-Cam observations (PI: D. J. Sand). 

2.2.3. Spectroscopic coverage 

MS2137 was targeted in two complementary 
spectroscopic campaigns using the Magellan Tele- 
scopes: 115 24 ^in-selected galaxies in the clus- 
ter core, observed using LDSS3 (PI: E. Egami); 
IMACS masks targeting Herschel sources (PI: T. 
Rawle), exploiting the larger field of view to ob- 
serve SPIRE detections beyond the LDSS3 field. 
The merged catalog contains 231 sources within 
the SPIRE footprint, of which 78 are in the clus- 
ter (0.30 < z < 0.33). Nine sources were observed 
by both instruments, with no disagreement. 



2.2. MS2137.3 2353 

MS2137.3-2353 (RA = 21:40:15.10, Dec = - 
23:39:39.0; z = 0.313; hereafter MS2137) is an 
ideal counterpoint to the Bullet Cluster. Located 
at similar redshift, the massi ve cluster is an ap- 
proxi mately spherical system ( Sand. Treu fc Ellis] 
2002), with an isolated BCG coincident with the 
X-ray center (a reliab le indicator of an undis- 
turbe d cluster; 
2009h . 



e.g. Sanderson. Edge fe Smith! 



2.2.1. Herschel observations 

MS2137 has HLS FIR coverage comparable to 
the Bullet Cluster SDP, with 8'x8' PACS maps at 
100 and 160 (xm and three-band SPIRE coverage 
(250-500 jum). The PACS sensitivities (3 a limits 
of 2.5 and 4.6 mJy) are marginally lower than in 
the Bullet Cluster SDP maps as the efficiency of 
PACS improved during the first year of operation. 
The nominal detection limits of the SPIRE maps 
are consistent with the Bullet Cluster (8.0, 9.3, 
10.3 mJy). All Herschel instrument configurations 
and data reduction are as described for the Bullet 
Cluster in Section |2~T1 



2.2.2. Additional imaging 

The spatial coverage in the Spitzer bands is 
more restricted for MS2137 than the Bullet Clus- 
ter, as shown by Figure [TJ However, Spitzer ob- 
servations did occur during the cryogenic phase 
(Pis: G. Rieke, E. Egami), so all four IRAC chan- 
nels and MIPS 24 /im are available, with the same 



2.3. Field samples 

We compile low-density comparison samples 
with the same observed photometry as the clus- 
ter population. These field samples comprise all 
galaxies spectroscopically confirmed to be in the 
foreground of either cluster (0.05 < z < 0.26). The 
most local galaxies (z < 0.05) are excluded as our 
filtered PACS maps are not optimized for extended 
sources. Background sources are also not included 
to avoid the additional uncertainty of lensing on 
the intrinsic luminosity. Although we mostly em- 
ploy these two foreground samples simultaneously, 
we maintain the distinction between the fields be- 
cause of the difference in photometric coverage 
(specifically a lack of 70 fxm for MS2137). The 
field and cluster sources are treated identically at 
every stage of analysis. 

W e also compare to large blank-field samples 
fro m iBlain. Barnard fc Chapman! (|2003l BBC03) 
and iHwang et al.l (|2010L H10). selecting the 129 
sources in the redshift range 0.05 < z < 0.5. 
BBC03 is a compilation of earlier studies using 
mostly IRAS data in the infrared. H10 com- 
prises two subsamples, an SDSS-selected sample 
observed by AKARI (predominantly z < 0.1), and 
a ffersc/id-detected sample in the GOODS-North 
field (mostly higher redshift). Both contribute 
sources in the required redshift range for this com- 
parison. 



See Sect ions [3 . 3 . 1 1 and 14 . 1 1 for more details 
on the suitability and usage of the external field 
samples. 
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3. Analysis 

3.1. Mult i- wavelength counterparts 

The final samples are derived by matching 
MIPS 24 /im and (optical) spectroscopic cata- 
logs. A 24 /zm detection (beam size 6") could 
be the counterpart to any optical imaging source 
within the MIPS rms pointing error (1.4"), and 
all were initially viewed as potential matches. The 
deep IRAC maps contain many more sources than 
MIPS 24 /zm, including a significant fraction of 
optical detections. The most likely optical coun- 
terpart to any MIPS source can be selected by in- 
teractive examination of flux and position through 
all available intermediate (i.e. near-infrared and 
IRAC) bands. The optical imaging counterparts 
are then cross-referenced to the spectroscopic cat- 
alogs to derive the matched MlPS-spectroscopic 
catalog. A conservative approach is adopted to en- 
sure MIPS sources are not assigned incorrect red- 
shifts; spectroscopically-confirmed optical sources 
with unsure matches in the IRAC or MIPS bands 
are not included in the final samples. 

Given the detection limits of infrared imaging 
available for the two clusters, typical SEDs show 
that any source at z~0.3 with a significant Her- 
schel detection, will also be detected at 24um. 
Identifying the corresponding sources in the un- 
confused PACS bands (beam size ~8") is straight- 
forward, and the small number of targets in this 
study allows for individual examination of every 
source. An automated or statistical counterpart 
algorithm may miss sources with unusual SED 
shapes. In the confusion-limited SPIRE bands, 



Table 1: Final sample sizes 



Total a Herschel b Low L 



IR 



Bullet Cluster 


64 


37 


37 


MS2137 


20 d 


11 


9 


Bullet field 


34 


20 


19 


MS2137 field 


15 


11 


9 



"Spectroscopically-confirmed cluster members with 24 (im 

detection; see footnote d 
'...and at least two Herschel detections for FIR SED fitting 
c ...and L IR < 2xlO n L (for Tabled 
"Includes three isolated spectroscopic cluster members 

detected by SPIRE beyond the MIPS 24 (im coverage 

(see Section 3.1) 



counterparts are not located directly, as flux is 
assigned to all MIPS sources simultaneously (de- 
scribed further in the following section). Three 
bright Herschel sources in MS2137, which lie be- 
yond the MIPS 24 /zm coverage are matched un- 
ambiguously to isolated optical sources identified 
in Subaru imaging. These sources are included 
despite the lack of 24 /im photometry; removing 
them from the sample does not affect the results. 

The final sample sizes are presented in Table 
[TJ with the cluster members also plotted in Fig- 
ure [H The Bullet Cluster analysis presented in 



iRawle et al.l ( 2010l ) comprised 23 sources, 14 of 
which had only SPIRE photometry. The new 
PACS imaging completes 70-160 /im coverage for 
all these sources, while also increasing the FIR- 
detected cluster sample by a further >50%. 

3.2. Photometry 

For Spitzer and PACS, where blending and 
source confusion is minimal, fluxes are mea- 
sured using simple ape rture photometry via th e 
Sextractor package ()Bertin fc Arnoutd Il996h . 
Aperture corrections are applied to each band, 
based on instrument handbook values. As rec- 
ommended for PACS imaging produced by the 
pipeline (and calibration product) used in this 
studjH, photometry is calculated in 12", 12", 22" 
radius apertures (at 70, 100, 160 /jm respectively) 
with aperture corrections of 0.886, 0.886, 0.916. 
Background subtraction is not required prior to 
aperture photometry in PACS maps, as the high- 
pass filtering artificially sets the mean background 
to zero. Flux uncertainties are estimated by Sex- 
tractor from the global sky rms error. 

In the SPIRE bands, where background source 
confusion is above the instrument noise, PSF fit- 
ting is the favored technique for measuring fluxes. 
We employ the Iraf routine Allstar, simulta- 
neously fitting a model PSF to all MIPS 24 /im 
positions. The use of a deep prior catalog removes 
a large fraction of the contaminating background, 
minimizing the need to de-boost SPIRE fluxes. 
Significant residual flux (e.g. from faint high-z 
sources below the 24 /im detection limit) would 
be identified by unphysical discontinuities in the 
final SED. In practice, for PACS sources at z<0.5 



2 PACS observers' manual; 

\http://herschel.esac.esaAnt/Docs/PACS/html/pacs.om.html\ 
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the detection rate in the SPIRE 250 /zm band is 
~50%. At 350 and 500 /im the Rayleigh- Jeans 
tail becomes very faint, and the detection rate is 
even lower at 20% and 5% (3 sources) respectively. 

3.3. Infrared SEDs 

3.3.1. SED fitting 

Infrared luminosity is derived via the best- 
fitting template to all MIPS and Herschel pho- 
tometry (i.e. A bs > 24 (Um). IRAC data is not 
included in the fits as these wavelengths are dom- 
inated by poorly constrained PAH features, and 
only sources with at least two Hersch e l dete ctions 
are analyzed. We us e thelRieke et all (|2009l R09) 
and ICharv fc Elbaz] ([200l|, CE01) templates, de- 
rived from a small number of star-formation dom- 
inated local (U)LIRGs interpolated to produce a 
library of galaxy spectra classified by IR luminos- 
ity. Template shape varies significantly with lu- 
minosity, in the sense that higher luminosity dust 
peaks at shorter wavelengths. As we are interested 
in studying this shape for cluster galaxies, we ig- 
nore the nominal luminosity class of the templates 
to find the best fit via a least squares minimiza- 
tion (from the Python SciPY. Optimize pack- 
age) , and allow a free-floating normalization factor 
to scale the template to the intensity of the source. 
Hereafter the luminosity class of the best-fitting 
template will be referred to as the "template lu- 
minosity" , ^template- 

The true total infrared luminosity of each 
galaxy (£ir) is derived by integrating the scaled, 
best-fitting template over the (rest frame) wave- 
length range Ao = 8-1000 /im. Uncertainty in the 
luminosity is estimated by Monte Carlo simula- 
tions of the fitting procedure. The star formation 
rate (SFRir) is calculated f rom Ltr fo ll owin g the 
simple relation presented in iKennicuttl (jl998l n. 

The characteristic temperature of the dust 
component is calculated by fitting a single- 
temperature greybody of the form 



is observed in the ra nge ~l-2 (e.g. = 1.8 in the 



S u = N^B v {T dust ) 



(1) 



where is flux density, A is a free-floating nor- 
malization factor, is dust emissivity index and 
B v (Tdust) is the Planck blackbody radiation func- 
tion for a source at temperature Tdust- Typically, 



solar neighborhood, iPlanck Collaboration et al 
20111 ). We fix = 1.5 to mirror the analysis of 



our external comparison samples BBCO3+HK0, 
although we note that — 2 is another popular 
choice in the literature. /3Td us t remains approx- 
imately constant, so using 0=2.0 would system- 
atically reduce the derived temperatures by 10%. 
Tdust for a z~0.3 cluster galaxy is usually well 
constrained, as the SED peaks between 100 and 
160 fim for dust at ~30 K. In contrast to the 
template fits described above, 24 ^m photometry 
is not included in the greybody fit. Sources with 
fewer than two Herschel detections are excluded 
from the greybody fit procedure. 

3.3.2. Robustness of SED fits 

First we compare the derived infrared luminos- 
ity from the best fit CE01 and R09 templates. The 
mean ratio L m (CE01)/L m (R09) = 0.90 ± 0.02, 
with an rms scatter of 0.08 dex. Lir is generally 
independent of template set, with a marginal ex- 
cess in the mid-infrared of the best fit CE01 tem- 
plate for sub-LIRG galaxies. Appendix IA1 expands 
this analysis by comparing Lir(R09) to values ex- 
trapolated from 24 /zm data alone. 

We now examine the dependence of our fits 
(and hence the derived values of the primary pa- 
rameters Lir and Td us t) on the availability of data 
at various wavelengths. We test the effect of non- 
detection in SPIRE bands by temporarily exclud- 
ing photometry at A > 250 /im. Lir and Tdust typ- 
ically differ by <0.1 dex and <1 K (respectively) 
despite the poorer constraint on the Rayleigh- 
Jeans tail. The difference is larger (~0.3-0.4 dex 
and 3 K) for sources with particularly cool charac- 
teristic dust temperatures (Td us t < 25 K), where 
the dust component peak lies closer to the SPIRE 
bands. Typically, the dust SED maximum lies 
within the PACS bands for sources at the clus- 
ter redshift and below, and therefore PACS data 
is the primary constraint on the best fit. 

A major difference between the Bullet Cluster 
and MS2137 data is the lack of 70 /xm coverage in 
the latter. For sources in the Bullet Cluster (and 
foreground field), we compare parameters derived 
from fits excluding 70 /jm photometry, to those 



3 SFR IR [M yr" 



4.5x10- 



-IR [erg s 



4 BBC03 and H10 employ exactly the same greybody formu- 
lation as we use, with /3=1.5 
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derived in Section 1331 LiR(no 70 /im)/ZiR = 1.0, 
with an rms scatter of 0.08 dex (for either tem- 
plate set). Figure [2] compares Tdust derived when 
including/excluding 70 /im photometry, and shows 
an increased scatter for sources with Td us t > 30 K. 
Although Si6o/<Sioo is a reasonable tracer of Td us t 
(see Appendix IBI for more details, especially Fig- 
ure [12]), above ~30 K, 70 /im photometry adds 
significant additional information about the SED 
shape. Excluding sources Td us t > 30 K from the 
previous Lin test reduces the rms scatter to 0.04 
dex (from 0.08 dex). These results suggest that 
70 /im data is required to fully constrain the peak 
of dust component in the warmest cases; without 
70 /im photometry, the uncertainty on Lir and 
Tdust in warm dust components is approximately 
doubled. We return to this bias in Section |4~2] 

3.3.3. Observed SEDs 

Figure [3] displays IR SEDs for a selection of 
sources in the Bullet Cluster. In all cases the best 
fitting template is a good description of the overall 
shape of the observed SED. The galaxies shown 
in the upper row are best fit by templates with 
a luminosity class consistent with the integrated 
luminosity (i.e. L temp i ate ~L m ). In other words, 
the shape of the dust component is the same as 
a local galaxy of similar luminosity. Although the 
best fitting template for the SEDs in the lower 
two rows generally correspond well with the ob- 
servations (the lower-left panel; marked WD4is an 
exception) , ^template is more than an order of mag- 
nitude higher than the integrated Lir. The shape 
of the dust component resembles a much more 
luminous local template. Figure 2] displays the 
template mi s match for all the Herschel sources. 
Rawle et al. I <|2010h attempted to explain the dis- 
crepancy via an excess in Sioo/S24- Appendix [B] 
presents a detailed exploration of this discrepancy. 
We suggest that the peak wavelength (traced by 
e.g. S160/S100) is a much more successful probe 
of the unusual behavior. Furthermore, Figure [4] 
shows that for both template sets, the discrepancy 
[i.e. log(L t em P iate/£iR.) ^ 0] is well correlated with 
characteristic dust temperature, a parameter the 
templates effectively keep constant at a given lu- 
minosity. Sources with the largest template dis- 
crepancy have unusually warm or cold character- 
istic dust temperatures for their luminosity. 
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Fig. 2. — The effect of 70 jiim data on the charac- 
teristic dust temperature Td us t for sources in the 
Bullet Cluster (filled red circles) and foreground 
field (open red circles). There is significantly more 
scatter for warmer sources, where 70 [im. is re- 
quired to constrain the peak. When 70 /im pho- 
tometry is excluded, sources below the dotted line 
no longer have sufficient data for a greybody fit, 
and are plotted at 14 K for convenience. 
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Fig. 3. — IR SEDs of spectroscopically-confirmed member galaxies in the Bullet Cluster. Uppermost row: 
Three examples of normal star-forming galaxies, for which Ltempiatc ~ -^iRi Lower rows: The six 'warm 
dust' galaxies with -^template > -^ir- In each panel, observed photometry is plotted in black. Best fitting 
templates are shown as solid lines (red=R09, blue=CE01), with LrR and ^template given in matching colors 
in the upper-left. The template with -k t cmpiatc=£iR ^ s added as a dashed line normalized to the same peak 
flux as the solid line. The best-fitting greybody is displayed as a grey dashed line, with the characteristic 
dust temperature shown to the lower-right. 
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Fig. 4. — The template mismatch, log(Ltempiatc/£iR.), compared to characteristic dust temperature (Tdust) 
derived from the best-fitting greybody (left panel=R09; right panel=CE01). Bullet Cluster=red filled circles, 
MS2137=bluc filled squares. Field samples in corresponding open symbols. Solid lines show the best linear 
fit to the data, implying a strong correlation between the parameters. The (yellow hatched) shaded regions 
guide the eye to sources for which |log(Ttempiate/£iR,)| < 0.5 (i.e. those without a significant template 
mismatch) . 



4. Results and discussion 

Throughout the remainder of this study, we 
use the infrared luminosity derived from R09 tem- 
plates. Our conclusions would not change if we 
instead adopted the CE01 templates, as noted in 
Section [5X2] 

4.1. IR luminosity temperature relation 

A relationship between dust temperature and 
infrared luminosity is often assumed, usually 
implicitly through the use of templates (e.g. 
Dale fc Heloul 20021 and t h ose in this study). 
Blain. Barnard fc Chaprnanl (|2003l BBC03) pre- 
sented the observed luminosity-temperature re- 
lation, concluding that the scatter was too large 
to estimate temperature uniquely from luminos- 
ity (or vice versa). BBC03 also highlighted a 
stark dichotomy in this parameter space. Whilest 
the brightest LIRGs (L m > 2xlO n L Q ) at all 
redshifts were observed to have Td us t in the range 
30-70 K, less luminous objects rarely have Td us t > 
40 K with a mean ^30 K. Consequently, the sub- 
LIRG population should be treated separately, 



and we adopt a cut at the turning point, Lm — 
2xl0 n L Q . 

We explore the distribution of our Herschel 
sources in the £iR-Td us t plane in Figure [5l Our 
field samples (open colored circles and squares) 
lie within the locus of 0.05 < z < 0.5 sources 
from BBC03 and H10 (small grey symbols). Two 
MS2137 field galaxies have Td us t~37 K, but also 
iiR > 4xl0 12 L Q , so are not unusually warm 
for their luminosity. One possible outlier is the 
Bullet field source with Td us t = 37 K and Lir ~ 
1.2xl0 10 Lq (HLSJ065807-555541). This galaxy 
falls within a minor foreground over-density at z 
= 0.234, so may be a group galaxy rather than an 
isolated field source. All other field galaxies in our 
samples have Td us t < 33 K. 

MS2137 cluster members (filled blue squares in 
FigureEJ) also lie within the BBC03+H10 distribu- 
tion. The two warmest sources in MS2137 (Td us t = 
33-37 K) are the most luminous cluster galaxies in 
this study, Lir, — 2-4xlO n L Q . Both are on the 
periphery of the cluster and neither are covered 
by PACS, so the dust peak is not well constrained 
(Tdust is likely to be a lower limit). No sub-LIRG 
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Fig. 5. — Main panel: IR luminosity - dust temperature relation for the cluster (red filled circles=Bullet, 
blue filled squares=MS2137) and field samples (corresponding open symbols). Lir is derived from the 
R09 templates, and Td us t f r om a single greyb ody. Grey symbols display 0.05 < z < 0.5 sources from 
Blain. Barnard fc Chaprnanl (|200l BBC03) and iHwang et all (|2010l H10). The vertical dashed line guides 
the eye to the dichotomy at Lir ~ 2xlO n L Q . Several galaxies in the Bullet Cluster have significantly 
warmer dust (Td us t > 37 K; above the horizontal dotted line) than similar luminosity field galaxies. All but 
one field galaxies with Lir < 2xlO n L© have Td us t < 33 K. Right panel: Normalized distribution of dust 
temperatures for Lir < 2xlO n L Q galaxies in the Bullet Cluster (red outline), Bullet foreground sample 
(magenta hatching) and BBC03+H10 field (grey solid). 
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sources in MS2137 have T dust > 33 K. 

In contrast, the Bullet Cluster (filled red circles 
in Figure[5]) contains several obvious outliers in the 
luminosity-temperature plane, visible as a high- 
temperature tail in the normalized histogram of 
sources with Lir = 2xlO n L Q , displayed in the 
right-hand panel of Figure [5] A similar tail is not 
identified in the corresponding field distributions. 

4.2. Significance of the warm dust popula- 
tion 

The number of Lir < 2xlO n L Q sources is 
small: 37 in the Bullet Cluster, nine in MS2137 
and 28 in the two foreground fields (19 of which 
are in front of the Bullet Cluster). We test whether 
the lack of galaxies exhibiting warm dust tempera- 
tures in MS2137 and the field samples is significant 
by calculating the likelihood of randomly drawing 
the observed sample from an identical tempera- 
ture distribution to the Bullet Cluster members. 
We confine this analysis to sources with Lir < 
2xlO n Lq, where the luminosity-temperature re- 
lation is flat. Table [2] summarizes the results, con- 
centrating on the probability of selecting the ob- 
served number of sources above three temperature 
thresholds (T dust = 33, 37, 40 K). 

The lack of warm dust hosts in MS2137 is not 
significant. As described in Sections [2] and El there 
are multiple inconsistencies between the datasets 
(mid-infrared footprint, 70 ^m coverage, depth at 
160 /im). Particularly important is the lack of 70 
jttm coverage (as discussed in Section l5.3.2[) , as the 
peak wavelength of a warm dust component falls 
below 100 /im. Warm dust sources in MS2137 may 
simply be misclassified without 70 /im photome- 
try to constrain the peak. Comparing the field 
samples to the Bullet Cluster, Table [5] reveals a 
marginally significant difference in the distribu- 
tions (> 3 a). For L m < 2xlO n L , the Bullet 
Cluster has more Td us t > 33 K sources than our 
foreground fields, and more extreme dust (Tdust 
> 40 K) than the BBC03+H10 sample (also il- 
lustrated by the histogram in Figure [5]) . Warm 
sources are at least preferentially located in the 
cluster environment. 

In summary, the dust temperature distribu- 
tion of the Bullet Cluster sources is not consis- 
tent with the foreground field samples. Here, we 
formally define 'warm dust' sources as those satis- 



fying the following two criteria: (1) Lir < 2xlO n 
Lq, marking the known break in the IR luminos- 
ity - temperature relation; (2) Tdust > 37 K, the 
nearest integer temperature to the 3 a tail of the 
BBC03+H10 temperature distribution. There are 
six warm dust sources in the Bullet Cluster: their 
SEDs are displayed in the lower two rows of Figure 
and derived properties are also listed in Table 
As dust temperature and redshift are degen- 
erate, it is worth noting that the six warm dust 
sources all have concurring spectroscopic redshifts 
from Magellan/IMACS and CTIO/Hydra, ruling 
out the possibility of foreground interlopers. 

4.3. AGN contribution 



We now explore whether obscured AGN are 
contributing to the mid-infrared of the warm dust 
sources, elevating the characteristic dust temper- 
ature. We note that AGN could only account for 
the variation between field and cluster samples 
if star-forming AGN hosts are more prevalent in 
clusters. At the redshift of the Bullet Cluster, clus- 
ters are likely to contain a lower fraction of AGNs 
than comparable field samples (a few per cluster; 
Martini et al1l2006t iGalametz et al.ll2009h . 



We examine several diagnostics to assess AGN 
contribution to the infrared. First we consult the 
Chandra point source catalog for the Bullet Clus- 
ter (private communication, M. Markevitch) , con- 
taining 145 sources in the central ~20' to a flux of 
2.5 xlO -16 erg cm -2 s -1 . Six Herschel-detected, 
spectroscopic sources are coincident (Ar < 1") 
with Chandra point sources; three in the fore- 
ground and three cluster members. Two of the 
cluster members have {Lir =2.7x10 10 L q , Td us t 
= 27 K} and {L m =2.4xl0 10 L , T dust = 30 
K} (L x = 4xl0 42 erg s _1 , 2xl0 41 erg s" 1 re- 
spectively). The third cluster member is a warm 
dust sources (WD1; f(0.5-2 keV) = 2.83xl0" 15 
erg cm -2 s _1 ; Lx = 8xl0 41 erg s -1 ). Based on 
X-ray non-detection, the remaining five warm dust 
sources have X-ray luminosities Lx < 7xl0 40 erg 
s -1 . Stacking on the optical positions of these 
five sources in the archival Chandra X-ray image 
reveals no significant signal, further constraining 
the average X-ray luminosity of the other warm 
dust hosts to L x < 4xl0 40 erg s _1 . 

A simple power-law continuum through the 
IRA C photometry (a sim ple indicator of an AGN; 
e.g. Donley et al. 20081) is not a good fit to the 
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Table 2 

Significance of the number of warm dust galaxies 



Lm < 2xlO n L Q 






ast>40 K 


Tdust 


>37 K 




as t>33 K 


low luminosity subsamples 


"sample 


n 


P 


n 


P 


n 


P 


Bullet Cluster 


37 


3 




6 




10 




MS2137 


9 





38% 





13% 





3.8% 


BBC03 & H10 field (0.05<z<0.5) 


46 


1 


<0.1% 


4 


1.7% 


9 


3.2% 


Bullet & MS2137 field 


28 





1.3% 


1 


0.1% 


1 


<0.1% 


Bullet field 


19 





10% 


1 


7.1% 


1 


0.6% 


Bullet field (excl. HLSJ065807-555541) 


18 





13% 





1.2% 





<0.1% 



Probabilities reflect the likelihood that the sample is drawn from an identical parent distribution to the 

Bullet Cluster. 



Table 3 

Derived properties of the warm dust sources in the Bullet Cluster 





full ID 


z 


L m (R09) 
[xlO 10 L ] 


SFR (R09) 
[M e yr" 1 ] 


Tdust 

[K] 


M^ust 

[xlO 8 M Q ] 


AGN? 


WD1 


HLSJ065817- 


-560249 


0.300 


5.5 ± 0.3 


9.5 ± 0.5 


39 ± 1 


2.1 ± 0.4 


(1) 


WD2 


HLSJ065826- 


-555913 


0.301 


1.9 ± 0.2 


3.3 ± 0.3 


43 ± 3 


4.4 ± 0.2 




WD3 


HLSJ065829- 


-555316 


0.304 


9.4 ± 0.7 


16 ± 1 


39 ± 2 


3.4 ± 0.6 




WD4 


HLSJ065831- 


-560336 


0.290 


6.8 ± 0.4 


12 ± 1 


48 ± 2 


1.0 ± 0.2 


(2) 


WD5 


HLSJ065842- 


-560041 


0.296 


2.5 ± 0.3 


4.3 ± 0.4 


37 ± 1 


1.2 ± 0.3 




WD6 


HLSJ065850- 


-560037 


0.289 


4.3 ± 0.3 


7.5 ± 0.5 


43 ± 2 


1.0 ± 0.2 





(1) coincident with X-ray point source; see Section l4~3l 
(2) AGN-like SV0/S24; right panel, Figure [6] 
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warm dust SEDs. In the left panel of Figure [6] we 
show an IRAC color-color dia gram for mor e detai l. 
Following the prescription of IStern et al 
normal galaxies (without AGN) occupy the region 
below and to the left of the bound area. Fore- 
ground and cluster sources appear separated due 
to the peak of stellar emission shifting red-ward 
through the IRAC bands with increasing redshift. 
Dominant AGN are located within the bound area, 
and none of the warm dust sources are found here. 
Galaxies with possible sub-dominant AGN contri- 
bution are found to the right of the AGN locus. 
All of the warm dust sources are found in this 
region, which may indicate a low level contribu- 
tion to the mid-infrared, or may be a consequence 
of an unusual warm dust component. The warm 
dust source coincident with an X-ray point source 
(WD1) exhibits the reddest IRAC colors of the 
war m dust pop ulation, and is not consistent with 



the IStern et al 



AGN criteria. 
The right panel of Figure [6] combines MIPS 
and PACS mid-IR d ata as a pro x y for the 
IR AS colors used by IVeilleux et all (|2009l ) (as 
in Hatziminaoglou et alT 2010 ). Here, dominant 



AGN have S 70 /S 2 4 < 8. CE01 and R09 template 
tracks indicate approximately where normal star 
forming-galaxies are located, and five of the six 
warm dust sources lie in this region (WD1 has 
SVo/SW = 13.1). WD4 is located below the AGN 
diagnostic line with S70/S24 ~ 2. However, this 
source does not resemble the normal AGN popu- 
lation which have much redder 5i6o/>Siou colors. 
Even if WD4 has a significant AGN component, 
the source does not resemble the typical AGN host 
population in the mid-infrared. 

4.4. Properties of warm dust sources 

4--4-1- Spatial distribution 

In terms of simple projected cluster-centric ra- 
dius, the warm dust sources are similarly dis- 
tributed to the whole FIR-bright population, i.e. 
absent from the densest core and found preferen- 
tially towards the periphery of the cluster, as illus- 
trated by Figure [7] Although spectroscopic fiber 
placement may contribute, this is consistent with 
the location of Herschel cluster sources in clusters 
with more comple te redshift data, e.g. LoCuSS 
([Smith et al.ll2010h . Five of the warm dust sources 
are located to the south of the cluster core, but 



the velocity offset ranges too widely (Acz ~ 4000 
km s _1 ) to be associated with a single substruc- 
ture (such as an in- falling group or filament). The 
warm dust sources generally avoid the very cluster 
center (>2-3 arcmin). Although they appear to be 
located preferentially in the south, this is due to 
the PACS coverage, which extends to larger radii 
in this direction (see Figure [TJ . 

The hot, dense intraclustcr medium (ICM) 
could plausibly trigger the mechanism responsi- 
ble for warm dust in some cluster sources e.g. by 
actively heating dust, or preferentially stripping 
cooler dust from the outer regions of the galaxy. 
Figure [7] displays the diffuse X-ray emission in 
the Bullet Cluster core, revealing that four of the 
warm dust sources are located in the vicinity of 
the densest X-ray gas, in regions where the gas 
temperature is > 10 KeV (for reference, the X- 
ray temp erature in the core of MS2137 is T x ~ 
5-8 KeV; ICavagnolo et al.ll2009t ). Two warm dust 
sources are located well away from the brightest 
X-ray emission (Tx ~ 8 KeV) , indicating that the 
extreme X-ray gas is unlikely to play a role. If the 
dense ICM were responsible for the warm dust, we 
may expect the phenomenon to be more prevalent 
in cluster galaxies than the observed ~15% level 
(6/37). 

4-4-%- Morphology 

If harassment or major interaction were respon- 
sible for the warm dust, the sources may have 
visible signs of disturbance. On the other hand, 
an undisturbed morphology would tend to favor 
a slow process such as strangulation (starvation) 
by the dense ICM. IMACS imaging (seeing ^0.5 
arcsec FWHM) hints at an irregular or disturbed 
morphology for three of the six warm dust galax- 
ies, although the resolution is too poor for in- 
depth study of galaxies at z^0.3. Only one warm 
dust source (WD4) is covered by high-resolution 
HST imagin g in th e core of the Bullet Cluster. 
Chung et all (|2010h investigated the source in de- 
tail and show it to be a normal barred spiral with 
no sign of major interaction (their figure 7). The 
relatively undisturbed morphology suggests that 



3 Guided by 24 /im flux alone lChung et al] (2OJLO5) classify the 
galaxy as the only ULIRG candidate in the Bullet Cluster; 
Herschel photometry reveals the source to be a sub-LIRG 
with unusually warm dust and an inflated 24 (im flux com- 
pared to Z/IR. 
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Fig. 6. — Infrared tests for dominant AGN contribution. Bullet Cluster=filled red circles, MS2137=fillcd 
blue squares (foreground fields are in corresponding open symbols). Warm dust sources (Td us t > 37 K) are 
highlighted by black star s and labelled with their WD number for reference (see Table [3]). Left panel: IRAC 
color-color plot with the Stern et"al] d2005l) AGN c riteria marked by grey dashed lines. Right panel: Mid-IR 
color-color diagram with the IVeilleux et al. I J2009h AGN diagnostic shown by the grey dashed line: sources 
below have >50% probability of a dominant AGN, and include one warm dust source (WD4). R09 (red) 
and CE03 (blue) template tracks are shown to indicate the location of normal star-forming galaxies, with 
log(Lt cm piatc/-k©)={10,ll,12} indicated by small numerals. 



strangulatio n is the dominant p rocess in, at least, 



this galaxy (jChung et al.ll201 



ttpr 

l3). 



4- 4- 3. Dust and stellar mass 

We briefly compare the dust properties to 
stellar mass, as the mass of the warm dust 
host galaxy could be an important parameter in 
stripping or interactions. Stellar mass is deter- 
mined using the SED - fitting method described in 
Perez-Gonzalez et al.l ( 2008I ). Dust mass is esti- 



mated from the greybody temperature and (rest- 
frame) 500 /im flux (calculated from the best fit 
greybo d y in S ection I3.3I) . using the formulation of 
iDraind ((20030. Dust masses are estimated from 
500 /im to minimize the dependence on tempera- 
ture. We note that the 500 fim flux may underes- 
timate the total mass by missing the warmest dust 
(which contributes only a fraction of the emission 
at this wavelengt h). Furt hermore, dust mass es- 
timates from the iDrainel models are wavelength 
dependent; dust mass derived from the flux at a 



(D 2 f 500 )/(K ahs B x (T dust )), where K abs (500 fim) 



rest wavelength of 200 ^m is systematically lower 
by ^0.2 dex, although the distribution of dust 
masses remains the same. 

Figure [8] plots dust mass against stellar mass, 
showing that the warm dust hosts tend to be 
amongst the least massive of our observed star- 
forming galaxies (4/6 have A/* < 10 10 M Q ). They 
are also all towards the lower end of the dust mass 
distribution, suggesting that the high character- 
istic temperature is not due to an influx of ad- 
ditional warm dust, but rather the conversion or 
removal of cooler dust. 

4-4-4- Interpretation of dust component SED 

Resolved infrared studies of nearby galaxies 
have investigated the internal gradients of the 
dust properties. Although the sample sizes are 
small, the general consensus is that bulges of 
late-type galaxies harbor enhanced dust heating 



( Alton et all Il998t iMelo et al.l l2002j) . Recently, 



Engclbra cht et al. [ (120101) explored the bulges and 



0.95 cm g 



discs of 13 nearby galaxies with Herschel, and con- 
firmed the general trend of warmer centers com- 
pared to the outskirts. Stripping preferentially re- 
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Fig. 7. — Location of the Bullet Cluster mem- 
bers with respect to the diffuse X-ray component 
detected by Chandra (blue intensity map). Her- 
shel detections are shown as filled red circles, with 
the warm dust sources highlighted by stars and 
labelled with their WD number for reference. 



10 



10 



1 1 



10 



t 10 9 



10 c 



10' 



: & Bullet 




o Bullet foreground field 


^ T dust > 


37 K 




• 


- © 


V* Q 
















10 9 10 10 10 11 
K [M ] 



10 



12 



Fig. 8. — Dust mass (calculated from the grey- 
body fit parameters) versus stellar mass (from a 
full optical/NIR SED fit) for the Bullet Cluster 
(filled red circles) and foreground field (open red 
circles) sources detected by Herschel. The warm 
dust sources are highlighted by stars as in previ- 
ous plots, and tend towards the lower end of the 
observed mass distribution. 



moves outer material, so the inferred characteristic 
temperature of an unresolved source would be seen 
to rise, as the cooler dust is removed first. The 
truncation of gas di sks in cluster galaxi es has long 
been observed (e.g. ICavatte et alJll994j ). while the 
suspected, associated truncation of the dust disk 
ha s recently been confi rmed using Herschel data 
bv lCortese et all The relatively low dust 

mass estimates of the warm dust sources in our 
study appears to support such a scenario. 

We investigate this further by examining 
whether the SED shape of warm dust sources 
can be recovered via the removal of cold dust 
from a 'normal' star-forming cluster galaxy. For 
the warm dust population and, separately, the 
remaining Bullet Cluster members, a composite 
SED is derived by summing the individual SEDs 
after normalizing at 24 /an. At the redshift of 
the Bullet Cluster, 24 /im traces the embedded 
star formation, which is most likely to remain un- 
changed by a stripping process. Figure [9] shows 
the resulting greybody fits to the two compos- 
ite SEDs, which recovers the population average 
parameters Td us t = 41 K and Tdust = 28.5 K re- 
spectively. For a given 24 (im flux, the warm dust 
sources have a smaller infrared dust component 
(and therefore lower dust mass) than normal star- 
forming cluster galaxies (as already observed in 
Figure [8|). The composite SEDs also suggest that 
an infrared luminosity extrapolated from 24 jj,m. 
would generally overestimate the true TrR for the 
warm dust sources, a result confirmed by Figure 
[TOl in Appendix El 

Assuming that the average warm dust source 
can be obtained by removing dust from the av- 
erage 'normal' star-forming galaxy, the stripped 
material equals the difference between the two 
composite SEDs. Subtracting the 41 K greybody 
from the other composite SED, leaves a 24.5 K 
cold component (demonstrated in Figure EJ . Us- 
ing a typical 24 /im flux (median of our Herschel- 
detected sample, S24 = 0.3 mJy), the dust mass of 
this cold component would be Mdust.coid ~ 8xl0 7 
Mq. The warm dust sources have a mean dust 
mass of ^5xl0 7 -5xl0 8 Mq, indicating that in 
this scenario, ~30-50% of the progenitor dust 
mass was stripped. Adding this amount of cold 
dust back into each of the individual warm dust 
sources, gives initial Tdust ~ 28 K and Ltr, ~ 
5xl0 10 - 2xlO n L<7), which would sit comfort- 
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• Composite warm dust source (best fit T dust =41 .0 K) 


o Composite "normal" source 
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Fig. 9. — Composite SEDs normalized at 24 /Ltm 
for the Bullet Cluster warm dust sources (Td us t > 
37 K), shown by solid circles, and 'normal' star- 
forming cluster galaxies, as open circles. Note that 
the 24 fj,m points for the two SEDs are coincident, 
due to the normalization. The range within each 
population is indicated by the grey bars. The best 
fitting greybodies are shown by thin black lines 
(T dust = 41 K and 28.5 K respectively). The nor- 
mal star-forming composite SED is also fit by a 
two component model, each shown by thick red 
dotted lines (total as red solid line): a warm com- 
ponent identical to the best fit warm dust source 
greybody, and a colder dust component (derived 
to have Td us t = 24.5 K). The dust in the lower 
temperature component would need to be removed 
to transform a normal star-forming cluster galaxy 
into a warm dust source. 



ably amongst the normal star-forming galaxy pop- 
ulation within the cluster (see Figure [5] for com- 
parison). These progenitors would also have dust 
masses towards the upper envelope of the normal 
star- forming galaxies (M dus t ~ 8xl0 8 - 8xl0 9 
M Q ; see Figure [8]) . It is not surprising that the 
detected warm dust sources are the progeny of 
galaxies lying at the high luminosity end of the ob- 
served range; less IR-luminous galaxies undergo- 
ing the same stripping process would have enough 
dust removed to render them undetectable in our 
Herschel imaging (Ltr < 10 10 L Q ). Stacking on 
-ffersc/ie/-undetected 24 /zm cluster sources recov- 
ers a marginal PACS detection at 70 /zm alone. 
Without a deeper FIR stack, we cannot confirm 
the presence of a lower-LiR, warm dust popula- 
tion. 

The stripping scenario appears plausible from 
the viewpoint of the total dust mass budget, via- 
bility of the progenitors, and the preferential strip- 
ping of lower temperature gas. However, the ob- 
served central dust temperature of nearby field 
galaxies is not as warm (<30 K) as the charac- 
teristic dust in some cluster sources (>40 K). The 
high dust temperature observed in several clus- 
ter galaxies appears to require additional heating, 
which may also be a consequence of the stripping 
mechanism (few of the resolved local field galaxies 
will have experience such a process). As a galaxy 
interacts with the ICM, central star formation may 
be triggered, in turn re-heating the surrounding 
dust. Stellar masses indicate that warm dust is 
located in less massive hosts, which suggests that 
molecular clouds in low mass galaxies are shocked 
and compressed as they encounter the dense ICM, 
triggering star formation in a relatively small vol- 
ume, which may resu lt in unusually warm dust 
(jBekki fc Couchl I2003T ). Heating and removal of 
outer (cold) dust are not exclusive, with both plau- 
sible components of a stripping mechanism that 
acts to increase the characteristic dust tempera- 
ture of a spatially-unresolved source. 



5. Summary 

We explore far-infrared bright cluster members 
from two dense environments at z^0.3, covered 
by the " Herschel Lensing Survey" : the merging 
Bullet Cluster and the more relaxed MS2137. At 
this redshift, characteristic dust temperature Td us t 
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is well correlated with the PACS color <S'i6o/<S'ioo, 
which straddles the FIR component peak. Pre- 
vious studies have shown that sub-LIRG galaxies 
display dust temperatures ~30 K, while dust in 
more luminous sources is warmer, Td us t = 30-70 
K. 

• Several warm dust galaxies are discovered in 
the Bullet Cluster, having higher than ex- 
pected characteristic dust temperature given 
their IR luminosity (Td ust >37 K at Lju < 
10 11 L ). 

• No warm dust galaxies are found in MS2137, 
although this could be an observational bias 
due to the lack of 70 /im coverage or the 
smaller FIR ficld-of-view. 

• Two field samples - the foreground of the 
Bullet observations, with identical photo- 
metric data to the cluster, and a larger lit- 
erature sample (|Blain. Barnard fc Chapman 
20031 : iHwang et alj|2010l) - are not compati- 



ble with the Bullet Cluster population, sug- 
gesting that warm dust galaxies are, at least, 
preferentially located in over-dense environ- 
ments. 

We examine scenarios which could be responsi- 
ble for the unusual population in the Bullet Clus- 
ter. One source is coincident with an X-ray point 
source (L x = 8xl0 41 erg s _1 ), but the remainder 
show no X-ray signature or optical/near-infrared 
line emission indicative of an AGN. The SVo / 5*24 
color of a second warm dust source suggests a sub- 
stantial mid-infrared AGN contribution, although 
the SED does not resemble a typical AGN- host. 

Composite SEDs for a warm dust galaxy and a 
'normal' star-forming cluster galaxy are derived 
to explore a simple dust stripping scenario, in 
which dust is preferentially removed from the out- 
skirts of a galaxy where dust temperature is lower. 
The mass of the removed dust is plausible, requir- 
ing 30-50% stripped from a normal star-forming 
galaxy to leave an SED similar to the warm dust 
sources. The central regions of (spatially-resolved) 
nearby galaxies are cooler than the characteristic 
temperature of the warm dust sources, suggest- 
ing that heating may be needed in addition to 
cool dust removal, a possible consequence of the 
same stripping mechanism. The only warm dust 



source with high-resolution imaging for morpho- 
logical analysis shows no sign of disturbance, fa- 
voring a slow stripping mechanism such as stran- 
gulation by the ICM (rather than e.g. harass- 
ment). However, we find no correlation between 
the location of warm dust sources and the densest 
ICM, indicating that the process may be triggered 
by the relatively less dense environment in the 
cluster outskirts, and suggesting that such sources 
should be observable in other, less massive clus- 
ters. 

With analysis of the full Herschel cluster sam- 
ple available from the lensing surveys, we will be 
able to better constrain the fraction of warm dust 
sources in the dense environment, and understand 
any dependence on environment, particularly the 
local ICM density, and host properties such as stel- 
lar mass and morphology. 

This work is partially based on observations 
made with the Herschel Space Observatory, a Eu- 
ropean Space Agency Cornerstone Mission with 
significant participation by NASA. Support for 
this work was provided by NASA through an 
award issued by JPL/Caltech. We would also 
like to thank the HSC and NHSC consortia for 
support with data reduction. This work has 
made use of the private version of the Rainbow 
Cosmological Surveys Database, which is oper- 
ated by the Universidad Complutense de Madrid 
(UCM). We would also like to thank Maxim 
Markevitch for providing the Bullet Cluster X- 
ray point source catalog, and Jean-Gabriel Cuby 
for the VLT/HAWK-I images. 
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A. Lir extrapolated from 24 fim 



Prior to the launch of Herschel, total infrared luminosity was often extrapolated from 24 /im photometry 
alone, where the mid-infrared sensitivity of Spitzer peaked. Initial Herschel observations show that such 
24 /im-derived IR-luminosities are reasonably reliabl e for star-formatio n dominated galaxies at low redshift 
(average scatter of 0.15 dex for sources at z < 1.5; Elbaz et al. 2010l ). We verify this for cl uster sources 
by calculating SFR 2 4 via equation 14 from Rieke et al.l ( 2009 ) , and use the Kennicutt (|l998l ) formulcQ to 
convert back to IR luminosity. After removing sources suspected of a dominant AGN contribution (those 
below the line in the right-hand panel of Figure 0, which are unaccounted for in the templates, L IR (24) and 
LiYi(Herschel) correlate reasonably well for our galaxies (rms scatter of 0.23 dex, Figure HU1) . It is difficult 
to constrain the AGN contribution for sources without 70 /im photometry. Including only confirmed star 
formation-dominated sources (i.e. excluding those without 70 /jm data), the rms scatter decreases further 
to 0.17 dex. The warm dust sources discovered in this paper (highlighted by squares in Figure [T0|). all have 
an over-predicted Lir from 24 /im. This suggests a mechanism that modifies the cold dust component while 
leaving the mid-infrared, which traces embedded star formation, unchanged (Figure [5]). 



B. Exploring the luminosity discrepancy 

For several Bullet Cluster galaxies, we identified a discrepancy between the nominal luminosity class 
of the best-fitting local template and the integrated lu minosity of that template when normalized to the 
observations (i.e. Lir ^ L tcmp i atc ). m Rex et al. ( 2010f) we discovered a similar (but opposite) discrepancy 



for high redshift sources, such that the shape of the SED resembles a low luminosity loc al galaxies (perhaps 



because star formation in local galaxies is more concentrated; i.e. iRuiopakarn |201ll) . Large blank- field 



surveys have not reported such a mismatch (up to z ~ 1.5; e.g. lElbaz et al.l I2010T ) . although this may 
be because an automated routine would report a good statistical fit, regardless of the nominal template 
luminosity class. The cluster sample presented in this paper displays a large range in this template mis- 
match, with no obvious correlation between L te mpiatc and Lir (Figures l4l and fTT|) . Note that the apparently 
smaller scatter for the R09 templates than CE03 is merely a consequence of a more limited range in L te mpiate: 
9.75 < log(L t emplate) < 13-0 compared to 8.4 < log(L tomplato ) < 13.5. 

Initial analysis of far-infrared data from HLS and the Local Cluster Substructure Survey (LoCuSS) sug- 
gested a large fraction of cluster galaxy SEDs were not w ell fit by local templates of comparable luminosity 



(|Rawle et aLll2010HPereira et al.ll2010l : ISmith et al.ll2010l ). The templates predict S100/S24 ~ 20-30 for L m 



< 10 11 galaxies at z~0.3 (Figure [T2l left panel), whereas the cluster sources show a wide variety of colors 
(5 < S100/S24 < 60). Although a low S100/S24 color can be attributed to an AGN increasing the flux at 24 
jum, larger ratios are harder to explain. These preliminary studies asserted that 25-50% of the star-forming 
cluster population had larger than expected S100/S24 (dubbed '100 /im excess'). However, S100/S24 tracks 
predicted by R09 and CE01 move in opposite directions with increasing luminosity, illustrating our uncer- 
tainty in this color space. The difference results from an elevated dust continuum level for higher luminosity 
templates in CE01, coupled to the deeper silicate absorption feature for higher luminosity templates in R09. 
The colors from this current study show a similar scatter around the template tracks (25% have S100/S24 > 
30; Figure fT2|) but there is no general correlation of S100/S24 with the luminosity discrepancy (Figure [T3l left 
panels). The scatter in S100/S24 is more likely to reflect different levels of strong silicate absorption and/or 
PAH features in the 24 /im band. S100/S24 is a poor probe of the template discrepancy. 

In contrast, both template sets display the general trend that increasing luminosity decreases Sieo/Sioo 
(Figure [T2l left panel). This is simply due to the FIR peak shifting blue- wards through the 160 and 100 /im 
bands, for sources at z~0.3. The right panels of Figure [13] show that as the peak shifts further from the 
mean wavelength (corresponding to Sigo/Sioo^l^), the mismatch between L te mpiatc and Lir becomes more 
pronounced. Figure |4] shows that the template discrepancy in our sample is highly correlated with Td us t (as 



7 SFR IH , [M® yr- 1 ] = 4.5xl0" 44 L m [erg sr 
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Fig. 10. — Comparison between 24 /xm-extrapolated Lir(24) and the Herschel Lir (via R09 templates). Left 
panel: Sources in the Bullet Cluster observations (both field and cluster members) are in red, MS2137 are in 
blue. The total rms scatter is 0.30 dex. Stars identify dominant AGN hosts (from Figure [5J right panel), a 
component unaccounted for in the simple 24 /im extrapolation. Excluding these, the rms scatter is reduced 
to 0.23 dex. The significance of any AGN component in sources without 70 fim photometry is unknown (open 
circles) . The remaining sources can be considered star- formation dominated (filled circles) . Right panel: The 
star-formation dominated sources replotted for clarity, with Bullet Cluster=red filled circles, MS2137=blue 
filled squares and field samples in corresponding open symbols. The rms scatter in this population is 0.17 
dex. The warm dust sources all have over-estimated Lir (24). Note that WD4 is the AGN host at the very 
top of the left panel, so does not appear in the right panel. 
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Fig. 11. — ^template compared to the integrated Lir of the normalized best- fitting template (see text for 
details): left=R09, right=CE01. Bullet Cluster=red filled circles, MS2137=blue filled squares. Field samples 
in corresponding open symbols. Note that the smaller apparent scatter for the R09 set is artificial as no 
templates exist log(L tC mpiatc) < 9.75. 
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derived from a greybody). S160/S100 is a useful tracer of Td us t for galaxies at this redshift (Figure [T2l right 
panel), and the correlation with template discrepancy is a consequence of the dependency on Tdust- 
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Fig. 12. — PACS color S160/S100 plotted against S100/S24 (left panel) and Td us t {right panel). Symbols and 
colors as in the previous plot. In the left panel, tracks of increasing luminosity are shown for the R09 (red) 
and CE01 (blue) templates at a redshift of z=0.3. The track positions for log(i te mpiatc)={10,ll,12} are 
indicated by small numerals. The right panel shows that S160/S100 is a reasonable probe of Td us t at this 
redshift (rms dispersion in temperature of 3.7 K). 
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Fig. 13. — The template mismatch, log(Ltempiate/im) (upper row=R09; lower row=CE01) compared to 
S100/S24 (left- panels) and S\eo/Sioo (right panels). The analogous plots for characteristic dust temperature 
(Tdust) are given in Figured] Symbols and colors as in the previous plots. The discrepancy between ^template 
and Lir is correlated with the FIR component peak, traced by <S'i6o/>S'ioo (at z~0.3) and driven in part by 
dust temperature. 
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